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Recently, emergent symmetry is one of fast-growing intriguing issues in many-body systems. Its roles and 
consequential physics have not been well understood in quantum phase transitions. Emergent symmetry of de¬ 
generate groundstates is discussed in possible connection to spontaneous symmetry breaking within the Landau 
theory. For a clear discussion, a quantum spin-1/2 plaquette chain system is shown to have rich emergent sym¬ 
metry phenomena in its groundstates. A covering symmetry group over all emergent symmetries responsible 
for degenerate groundstates in the plaquette chain system is found to correspond to a largest common sym¬ 
metry group of constituent Hamiltonians describing the plaquette system. Consequently, this result suggests 
that, as a guiding symmetry principle in quantum phase transitions, degenerate groundstates are induced by a 
spontaneous breaking of symmetries belonging to a largest common symmetry group of continent Hamiltonians 
describing a given system but can have more symmetries than the largest common symmetry. 

PACS numbers: 


I. INTRODUCTION 

Symmetry plays an important and indispensable role in un¬ 
derstanding physics. Revealing many properties of nature are 
allowed, in fact, owing to understanding a mechanism of sym¬ 
metry breaking M- As is known, symmetry of nature can 
be broken in two different ways. One is explicit symmetry 
breaking, e.g., the isotopic symmetry of the nuclear force j^]. 
The other is spontaneous symmetry breaking, which plays a 
more profound role than explicit symmetry breaking. As a 
modern terminology, spontaneous symmetry breaking, though 
appeared first in Baker and Glashow’s paper ||3], is consid¬ 
ered as an emergent phenomenon because once it happens to 
a system, an underlying physics of the system is not reducible 
to some sort of sum of behaviors of its parts, and not pre¬ 
dictable 0] • Spontaneous symmetry breaking phenomena 
have been observed ubiquitously in macroscopic systems in¬ 
cluding not only classical systems |0 but also quantum sys¬ 
tems such as quantum Ising systems, and also theoretically in¬ 
troduced quantum systems and experimentally prepared sys¬ 
tems that can be described by an effective Hamiltonian, e.g., 
spontaneous particle-hole symmetry breaking in the v = 5/2 
fractional quantum Hall effect Quantum phase transitions 
can then be a prototypical example for spontaneous symme¬ 
try breaking in many-body systems. The Landau theory for 
spontaneous symmetry breaking has guided how to under¬ 
stand quantum phase transitions 111 [j, S UOl ■ 

In the Landau theory, a consequence of spontaneous sym¬ 
metry breaking is a degeneracy of groundstates in broken- 
symmetry phase. Degenerate groundstates from spontaneous 
symmetry breaking is normally believed to have a less sym- 
me^ than their Hamiltonian for a given system parameter 
12,0 • More precisely, if a given system Hamiltonian H is in¬ 
variant under an unitary transformation U, i.e., UHU^ - H, 
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and the unitary transformation U is related to an element of 
the Hamiltonian symmetry group G, there are two possible 
groundstates, i.e., and U\tl/gs), that satisfies H\if/gs) = 
Egsl^ftgs) and HU\figs} = EgsU\figs), respectively. When the 
two groundstates become 4^ U\il/gs}, the system under¬ 
goes a spontaneous breaking of which symmetry consists of 
a subgroup g of the Hamiltonian symmetry group G in as¬ 
sociation with the unitary transformation U. Then the two 
degenerate groundstates seem to have a lower symmetry than 
the Hamiltonian. However, in general, groundstates can have 
a symmetry which does not belong to their Hamiltonian sym¬ 
metry. For instance, the two-fold degenerate dimerized states 
can have a U(l) symmetry induced by a local projector on 
the singlet state of the bond at the Majumdar-Ghosh point in 
the frustrated antiferromagnetic Heisenberg chain, so called, 
the J1-J2 model ifTTl] . Such a symmetry, which is absent in 
Hamiltonian for a fixed parameter, is called emergent sym¬ 
metry flll - i^ . However, the two dimerized states are not 
connected by a unitary transformation relevant to the emer¬ 
gent U(l) symmetry. The actual broken-symmetry giving rise 
to the dimer phase is the one-site translational symmetry be¬ 
cause the two dimerized states are not one-site translational 
invariant but one of them becomes the other under the one- 
site translational transformation. The broken-symmetry for 
the two degenerate groundstates belongs to the symmetry of 
the frustrated antiferromagnetic Heisenberg chain Hamilto¬ 
nian. In contrast to the normal belief, then, for a spontaneous 
symmetry breaking, induced degenerate groundstates do not 
have a broken-symmetry belonging to Hamiltonian symmetry 
but can have more symmetries not belonging to Hamiltonian 
symmetry. 

However, such an emergent symmetry phenomenon in 
groundstates can allow a more crucial question, i.e., are de¬ 
generate groundstates induced by the spontaneous symmetry 
breaking in the Landau theory? Supposed that there are two 
degenerate groundstates |i/ri) and Itf/i), i-S-, |tAi) 4^ in a 
system Hamiltonian H for a fixed system parameter, where 
the two groundstates satisfy = £’gj|i/^i/ 2 ). The two 



2 


groundstates can be related by a unitary transformation U, i.e., 
\4i\) - U \ip 2 )- The unitary transformation U can be applied to 
see how the system Hamiltonian can be transformed by U, i.e., 
UHU^. In general, due to occurring an emergent symmetry, 
there are two possible situations, i.e., (i) the Hamiltonian H is 
invariant under the unitary transformation, i.e., UHU^ - H, 
or (ii) the Hamiltonian H is not invariant under the unitary 
transformation, i.e., UHU'' + H. Straightforwardly, the case 
(i) is equivalent to what the Landau theory normally states 
about the spontaneous symmetry breaking. For the case (ii), 
however, the two degenerate groundstates might not be un¬ 
derstood within the spontaneous symmetry breaking in the 
Landau theory. One may then first ask whether such a situ¬ 
ation in (ii) can occur in any quantum system. In fact, the 
answer to this question is yes, i.e., the model system in Eq. 
([TJ will be shown to have such a two-fold degenerate ground- 
states, e.g., in the anti-ferromagnetic plaquette (AFP) phase 
{h + 0) and the staggered bond (SB) phase. Hence, it is nat¬ 
ural to ask whether an emergent symmetry can be introduced 
to understand why the degenerate groundstates occur in asso¬ 
ciation with the unitary transformation U. Further, can such 
degenerate groundstates be understood within a spontaneous 
symmetry breaking and Hamiltonian symmetry? 


For clear discussions on an emergent symmetry of ground- 
states in quantum phase transitions, we will numerically in¬ 
vestigate the spin-1/2 infinite plaquette chain in Eq. 

The infinite matrix product state (iMPS) representation yl 
1^ is employed for wavefunctions and the infinite time- 
evolving block decimation (iTEBD) method is used to 
get a groundstate wavefunction (See Fig. [TJ. To detect de¬ 
generate groundstates for a given parameter, we have used 
the quantum fidelity 12^ - 1^ defined as a overlap measure¬ 
ment between calculated groundstates and an arbitrary refer¬ 
ence state. From the singular behaviors of the quantum fi¬ 
delities, also, the phase boundaries are determined. To see 
symmetry of groundstates properly, we obtain an explicit form 
of groundstates from careful analysis on local magnetizations 
and two-site spin correlations in each phase. In order to under¬ 
stand a relation between a degenerate groundstate and Hamil¬ 
tonian symmetry, we investigate groundstates in whole pa¬ 
rameter range of the model Hamiltonian. Seven phases with 
two-fold degenerate groundstates and two phases with a single 
groundstate are clarified. The single groundstates in the two 
phases are found to have emergent symmetries. Furthermore, 
we find that a two-fold degenerate groundstates could not be 
understood within the spontaneous symmetry breaking pic¬ 
ture in Landau theory. Based on the common properties of the 
emergent symmetries in our model, we discuss and suggest a 
possible extension of the spontaneous symmetry braking pic¬ 
ture in association with a largest common symmetry group of 
the constituent Hamiltonians. Also, the spin structure factors 
for the degenerate groundstates with an emergent symmetry in 
each phase are shown to have experimentally distinguishable 
peak structures. 



FIG. 1: (color online) Top: Spin-1/2 plaquette chain in the quasi- 
one-dimensional lattice. The dotted box indicates a plaquette lattice 
where the four spins interact one another with the Ising interaction 
strength 7' denoted by the dashed lines. The thick solid lines indi¬ 
cate the Heisenberg interactions J between two spins. Bottom: One¬ 
dimensional spin chain mapped from the quasi-one-dimensional lat¬ 
tice for the infinite matrix product state (iMPS) representation. The 
labels from A-H' indicate a 16-site unit cell in the iMPS representa¬ 
tion. The labels in the one-dimensional lattice correspond to ones in 
the quasi-one-dimensional lattice. 


II. QUANTUM SPIN-I/2 PLAQUETTE CHAIN 

In order to investigate a relation between spontaneous sym¬ 
metry breakings and degenerate ground states, we consider a 
quantum spin-1/2 plaquette chain in a magnetic field in Fig. 
[1] The chain Hamiltonian with the Heisenberg intradimer and 
the Ising interdimer interactions can be written as 

^ — ^dimer + Hplaq Hfield 1 (1) 

where 

Hdimer ~ 72(S,-,„-S,-rf + S,-,-S,+i,,), (2a) 

i 

Hpia, - (2b) 

i 

Hf^eid = -hYfnu^sid+ni^nd- ( 2 ^^) 

i 

Here. Si^u,Si^d,Sij, and represent the spins in the /th pla¬ 
quette, where u, d, I, and r denote their corresponding lattice 
sites, i.e., respectively, at the up, down, left, and right sites of 
the plaquette (See Fig. m. Hdimer and Hpiaq describe the in¬ 
teractions on dimers and plaquettes with their corresponding 
interaction strengths J and 7' respectively. In Hfield, h is the 
strength of external magnetic field along the z direction. 

For J ^ J[ and h - 0, for instance, the Hamiltonian be¬ 
comes H X Hpiaq. Because Hpiaq consists of four-site Ising 
plaquettes with a periodic boundary condition and possesses 
a global Z 2 ® U(l) symmetry, the system can have degener¬ 
ate ground states due to a spontaneous Z 2 symmetry breaking 
in each plaquette. Also, for 7 » 7' and h - 0, the Hamil¬ 
tonian becomes the dimer Hamiltonian H » Hdimer that has 
a global SU(2)-rotational symmetry and then the dimerized 
spins are in a singlet state for 7 > 0. Thus the system can 
have a single ground state that has the same SU(2) symmetry 
with the dimer Hamiltonian. If the applied magnetic field is 
relatively very strong, i.e., /? » 7, 7', the system Hamiltonian 
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has a global U(l)-rotational symmetry and the system may 
be in a fully polarized ground state with the same U(l) sym¬ 
metry. Then, if the interaction parameters varies the system 
may undergo other types of spontaneous symmetry breakings 
because the above examples show that groundstates can have 
a different symmetry for different system parameters. In this 
study, for whole system parameter range, we will numerically 
investigate all degenerate ground states and their relations to a 
spontaneous symmetry breaking in the dimer-plaquette chain. 

For our numerical study, we use the IMPS algorithm ES. 
In order to employ the iMPS representation for wavefunc- 
tions, one needs to map the dimer-plaquette chain to a one¬ 
dimensional chain. In the mapped one-dimensional chain in 
Fig. [U the mapped Hamiltonian is four-site translational in¬ 
variant. To ensure a lattice symmetry breaking in the one¬ 
dimensional lattice, our numerical simulations are performed 
in the 4-site, 8-site, and 16-site iMPS wavefunctions. In fact, 
we have found that the 8-site and 16-site unit cells in the iMPS 
representations give a same result with a negligible numerical 
accuracy difference. Then, we will discuss our results based 
on the 16-site iMPS representation. 


III. FIDELITY PER LATTICE SITE AND DEGENERATE 
GROUND STATES 


For given parameters in an initial iMPS wavefunction, once 
the energy of the system becomes saturated during the iTEBD 
procedure, one can obtain a groundstate wavefunction. If a 
different choice of initial states leads to a different ground- 
state with a same saturated energy, the system has a degener¬ 
ate groundstate. When the system has more than one ground- 
states, degenerate groundstates can be distinguished by using 
the quantum fidelity per lattice site (FLS) with an arbitrary ref¬ 
erence state in Ref. |23. The FLS d{\ip), \(p)) ll^l^ is defined 
as 


lnt/(|i/r), |0» = hm ---, (3) 

L—»oo J_, 

where L is the system size. Here, the fidelity Fi\^), \(p)) be¬ 
tween a reference state \(p) and an iMPS ground state li/r) is 
defined as F{\ip), \(p)) — Kt/'|<^)|, where |i/f) e with 

n - 1, • • ■ ,N when the system has A-fold degenerate ground 
states As a scaling parameter ll^ in the thermody¬ 

namic limit, the FLS d is well defined with the characteris¬ 
tic properties, i.e., (i) normalization £/(|i/'), \(p) = \ip)) - 1 
with F(|i/r), 10) s 10)) = 1, (ii) symmetry ^f(|0), |0)) = 
^f(|0), 10)) with F(|0), 10)) = F(|0), 10)), and (iii) its range 
0 < d{\ijj), 10)) < 1 with F(|0), 10)) 6 {0,1). Actually, the FLS 
corresponds to a projection of degenerate groundstates into a 
chosen reference state. Then, a number of different values of 
FLS indicates the degeneracy of the groundstates. 

In our numerical calculation, for a given parameter, we 
have used many different initial states randomly chosen nu¬ 
merically to determine whether the system has a degenerate 
groundstate with the LLS. If only one value of LLS is detected, 
the system has a single groundstate, while if N different val¬ 
ues of LLS are detected, the system has A-fold degenerate 
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FIG. 2: (color online) Fidelity per lattice site for (a) h = 0, (b) 
h = 0.57 and (c) h = 1.57. In some ranges of 7'/7, FLS has one 
or two values. A single value of FLS indicates a single groundstate, 
while two different values of FLS indicate doubly degenerate ground- 
states. Also, the discontinuous FLSs at some values of 71/7 indicate 
occurring of phase transitions at those values. 


groundstates. However, if a chosen reference state is one of 
degenerate groundstates, the FLS cannot distinguish all of de¬ 
generate groundstates ll^ . If a chosen reference state is not 
one of degenerate groundstates, the FLS distinguishes all de¬ 
generate groundstates for a given system parameter. Thus, to 
distinguish all different degenerate groundstates properly, we 
have randomly chosen several reference states numerically in 
our FLS calculation. Also, we have calculated groundstates 
with different truncation dimensions = 4, 8, 16 and 32. 
The different truncation dimensions have been found to give 
a same numerical result. We present our numerical results for 
the truncation dimension= 32. 

In order to show clearly how to determine a phase boundary 
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from a characteristic behavior of FLS with degenerate ground- 
states, as examples, we display FLSs as a function of 7'/7 for 
(a) /i = 0, (b) /! = 0.5J and (c) h = 1.57 in Fig. El In 
Fig. EJa), we plot FLS as a function of 7'/7 for h - 0 with 
-37 < 7' < 37. For 7' < -1.4147, FLS has two values. 
At 7' = -1.4147, the FLSs jump into a single value of FLS. 
The observed single value of FLS for -1.4147 < 7' < 1.4147 
means the system has a single groundstate. Normally, a sin¬ 
gle groundstate implies that a spontaneous symmetry break¬ 
ing may not happen to the system. The further incensement 
7'/7 gives rise to two values of FLS for 7' > 1.4147. For 
7' < -1.4147 and 7' > 1.4147, the two values of FLS means 
a two-fold degenerate groundstates indicating that possibly a 
spontaneous symmetry breaking occurs. Note that the jumps 
of FLS occurs at 7' = -1.4147 and 7' = 1.4147. Such discon¬ 
tinuous behaviors of FLS indicate discontinues phase transi¬ 
tions ll^ E^l ■ In the parameter range, then, there are two 
phase transitions at 7' = -1.4147 and 7' = 1.4147. 

Similarly, in Fig. Etb), we plot FLS as a function of 7'/7 
for h - 0.57 with 0.57 < 7' < 27. As the interaction 7'/7 in¬ 
creases from 7' = 0.57, FLS has one value until 7' = 1.1187. 
The further incensement 7' /7 gives rise to two values of FLS. 
At 7' = 1.6327, the two values of FLS jump into another two 
values of FLS. Then, one can see that FLS has a single value 
for 0.57 < 7' < 1.1187 and two values for 7' > 1.1187. Note 
that at 7' = 1.6327 there occurs a phase transition between the 
two phases, both of which have two-fold degenerate ground- 
states. 

In Fig. EIc), we display FLS as a function of 7'/7 for 
h - 1.57 with 0 < 7' < 37. When 7'/7 varies from 0 to 
0.5, there is only one value of FLS corresponding to a single 
ground state. With a jump at 7' = 0.57, the FLS is split into 
two values for 7' > 0.57, which means a two-fold degener¬ 
ate groundstate. At 7' = 1.1187 and 7' = 2.597, as the 7'/7 
increases further, two values of FLS jump into another two 
values. Then, the three jumps at 7' = 0.57, 7' = 1.1187, and 
7' = 2.597 respectively indicate discontinuous phase transi¬ 
tions separating four different phases. 

Actually, by using the property of FLS, we have detected 
seven phases and their phase boundaries for the spin-1/2 
plaquette model. Figure |3] shows the phase boundaries be¬ 
tween the nine different phases denoted as singlet-dimerized 
(SD), fully polarized (FP), staggered bond (SB), modulated 
anti-ferromagnetic plaquette (MAFP), modulated ferromag¬ 
netic plaquette (MFP), staggered anti-ferromagnetic plaque¬ 
tte (SAFP), staggered ferromagnetic plaquette (SFP), anti¬ 
ferromagnetic plaquette (SAF), and ferromagnetic (F) phases. 
The SD and FP phases have a single groundstate from a single 
value of FLS. While the other seven phases, i.e., SB, MAFP, 
MFP, SAFP, SFP, AFP, and F phases, have a two-fold degen¬ 
erate groundstate. Note that, without knowing whether any 
spontaneous symmetry breaking occurs or not, our FLS is 
shown to determine the phase boundaries and the degeneracy 
of groundstates for each phase. For the parameter range, i.e., 
7,7', h > 0, a transfer matrix method has been used to ob¬ 
tain a similar phase diagram and groundstates in Ref. EH In 
order to discuss about a relation between spontaneous symme¬ 
try breakings and degenerate groundstates for the phases, one 
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FIG. 3: (color online) Groundstate phase diagrams for the spin-1/2 
plaquette chain for (a) 7 > 0 and (b) 7 < 0 in the J'.-h plane. The 
phase boundaries are determined by the discontinuous behaviors of 
the fidelity per lattice site (FLS). For each phase, its characteristic 
property is discussed in the text in association with relevant symme¬ 
tries of groundstates and the system Hamiltonian. 


should know an explicit symmetry of groundstates. To do this 
in our study, we will investigate a local property of ground- 
state wavefunctions, i.e., local mangetizations and two-point 
spin correlations for each phase in the following sections. 


IV. GROUNDSTATE WAVEFUNCTIONS AND EMERGENT 
SYMMETRIES 

Symmetries of calculated IMPS groundstates can be un¬ 
derstood from their local property and two-point spin corre¬ 
lations. Especially, local magnetizations can give us an in¬ 
formation about spin rotational symmetries of groundstates. 
Also, if they have a periodic structure, lattice symmetries of 
them could also be understood. Based on such properties of 
local magnetizations and two-point spin correlations, one may 
discuss symmetries of groundstate wavefunctions. However, 
the best way to discuss symmetries of groundstates wavefunc- 
tion requires knowing an explicit form of groundstate wave- 
functions. In our model, we have found that for our whole 
parameter range, i.e., all nine phases, only z-components of 
local magnetizations are nonzero from all groundstates, i.e., 
{S x) - 0 = (Sy). Furthermore, two-point spin correlations 
defined as C““'(|i - j\) - with a,a' e {x,y,z} are 
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FIG. 4: (color online) (a) Local magnetization (S^) at the lattice site 
i (i = 1,2, • • • correspond to A-H in regular sequence in Fig. □) 
and (b) two-point spin correlations as a function of lattice 

distance between two sites A and A-)- r for 7' = 0.27 and /z = 1.57. (c) 
Pictorial representation of the ground state with spin configuration in 
the FP phase. The circles in the diagram are spins at each site. A red 
arrow indicate a fully polarized magnetization. 


the spin states for the two sites are a product state. Also, if 
any pair of two sites satisfies the conditions in Eqs. ( l4al l and 
(Ebli, the state for the system is in a product state of the spin 
states of each site. For our FP groundstates, any two-point 
spin correlation satisfies the conditions in Fqs. (l4al l and (l4bl l. 
which implies that the groundstate is a product state of the 
spin states of each lattice site. Consequently, in the original 
quasi-one-dimensional lattice, with the fully polarized mag¬ 
netizations {Sp = {Sp - 1/2, the groundstate for the FP 
phase can be written as 

i 

In Fig. ID (c), the groundstate for the FP phase are presented 
pictorially. 

The groundstate in Fq. (|5]l has a U(l) rotational symmetry 
in the x-y plane. Also, the Hamiltonian in Fq. ([T]i has the U(l) 
rotational symmetry. Both the groundstate and the Hamilto¬ 
nian are one-plaquette translational invariant. Hence, for the 
FP phase, no spontaneous symmetry breaking occurs, which 
results in the single groundstate. As a result, this FP phase 
can be explained within the Landau’s spontaneous symme¬ 
try breaking theory. By comparing with the symmetry of the 
Hamiltonian, however, one can notice emergent symmetries 
for the groundstate. For instance, for a lattice rotation of each 
plaquette, the groundstate is invariant but the Hamiltonian is 
not. For a vertical-to-horizontal site exchange in each plaque¬ 
tte, e.g, (A, B) <-4 (C, D), also, the groundstate is invariant but 
the Hamiltonian is not. These imply that the lattice-rotation 
and the exchange symmetries are emergent for the ground- 
state. Consequently, the groundstate has such emergent sym¬ 
metries not belonging to the Hamiltonian symmetry. 


found to have a periodic structure. These facts allow us to ex¬ 
tract an explicit form of groundstate wavefunctions from the 
two-point spin correlations and local magnetizations for each 
phase. We will discuss emergent symmetries of groundstates 
in each phase. 

A. Fully polarized phase 

For the FP phase, there is a single groundstate from the FLS 
calculation. In Fig. St a), the local magnetizations at the lattice 
site i and (b) the two-spin correlations as a function of lattice 
distance r are plotted for 7' = 0.27 and h = 1.57. The local 
magnetizations have their maximum value, i.e. {S - 1/2 at 
all lattice sites. In Fig.|4jb), the two-point correlations satisfy 
{S ^-Sp — {S p{Sp for any pair of two spins. Actually, we have 

also observed that {SJSp ^ 0 ^ {S'.Sp and = 0 

with a, a' e {x,y,z} (not displayed). In general, if spin corre¬ 
lation between two sites i and j satisfies 

= <5f><5p, (4a) 

<5 = 0, (4b) 


B. Singlet dimerized phase 

Similar to the FP phase, for the SD phase, a single ground- 
state is detected from the FLS calculation. In Figs. |5ja)-(c), 
the local magnetizations at the lattice site i and the two-spin 
correlations as a function of lattice distance r are respectively 
plotted for 7' = 0.67 and h - 0.57. For the two-point spin 
correlations, we observe that {SfSp have a same behavior 
and {S“Sp*“) - 0. In Figs. |3b) andj^^c), thus, the two- 
point spin correlations (5 “5 pp show that except for the near¬ 
est two-spin correlation and all other two- 

point spin correlations are zero. In general, if two-point spin 
correlations between sites i and j satisfies 

^ {SfXSpp, (6a) 

for /^i+1, (6b) 

(S^Sp) = 0, (6c) 

the two spins for sites i and i + 1 are dimerized. Hence, the 
two sites (A, B) and (C, D) are dimerized. Also, the ground- 
state is a product state of the dimerized spin pairs, i.e., the two 
sites (A, B) and (C, D). Furthermore, the local magnetizations 
are zero at all lattice sites, which implies that the two bases 
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FIG. 5: (a) Local magnetization (S,) at the lattice site i for 7' = 0.67 
and h = 0.57 in the SD phase. Two-point spin correlations (b) 
<5^5 and (c) (5 g5 with a = x,y,zas a function of lattice dis¬ 
tance between two sites A/B and A/B + r for 7' = 0.67 and h = 0.57. 
For the singlet dimerized phase, two spins AB are dimerized and in 
a singlet pair state. Note that for the sites C and D, the two-point 
spin correlations <5 = (SA/a^A/B+r^- (d) Pictorial repre¬ 

sentation of the ground state with spin configuration. The circles in 
the diagram are spins at each site. A thick green line connecting two 
spins represents a singlet state. 


lit) > li)) at each site equally contribute for the local magneti¬ 
zations. The locally dimerized spin pairs then are in an oppo¬ 
site spin state, i.e., for instance, ITa) Ub) or |J,a) ITb) because 
= -1/4. This fact means that the two-spin state can 
be written by a linear combination of the two possible states, 
i.e., \tf/AB) = a ITa) \iB)+b IIa) ITb) with \a\^ + \b\^ = 1, where a 
and b are numerical coefficients. In this case, the zero magne¬ 
tization at each site give a constraint condition, i.e., |a| = \b\. 
Such a condition \a\ - \b\ allows only \a\ = |7>| = 1/ V2 with 


the normalization condition \a\^ + \b\^ = 1. In addition, from 
the property of {Sa ■ S^) = -3/4 with = -1/4, one 

can determine the coefficients explicitly, i.e., a — -b - \ j V2. 
Then, the dimerized spin pairs are in a spin singlet state (de¬ 
noted by a thick line between a pair of two spins in Fig.|5ld)), 
i.e., \\j/AB) = (ITa) Ub) - lU) ITb)) / V2- Similar discussions 
should lead to a spin singlet state for other pairs of two spins 
C and D. Consequently, for the SD phase, the groundstate can 
be written as 

l«A) - n 5 - |i,>T,+i,/>). 

(7) 

Now, let’s discuss a symmetry of the groundstate for the 
SD phase with h > 0. In the quasi-one dimensional plaquette 
lattice, the groundstate in Eq. © for the SD phase is one- 
plaquette translational invariant. Also, the groundstate in the 
SD phase has a SU(2) spin-rotational symmetry for h > 0, 
while for h - Q {J[ 0), due to the interdimer Ising interac¬ 

tions, the system Hamiltonian possesses a Z 2 spin-flip symme¬ 
try for the z-direction and a U(l)-rotational symmetry along 
the z-axis, and for h + 0, the original Z 2 ® U(l) symmetry 
of the Hamiltonian is broken into a U(l)-rotational symmetry 
due to the magnetic field applied to the z-direction. Thus, in¬ 
terestingly, in both cases h - 0 and h + 0, the groundstate 
wavefunction possesses more symmetries than the Hamilto¬ 
nian, i.e., Z 2 ® U(l) c SU(2) and U(l) c SU(2), respectively. 
Actually, the SU(2)-rotational symmetry of the groundstate 
could not be expected from the system Hamiltonian for the 
parameter range of the SD phase because the system Hamilto¬ 
nian is not invariant under the SU(2) rotation. In addition, the 
groundstate has a vertical-to-horizontal site-exchange sym¬ 
metry, i.e., the groundstate is invariant for exchanging sites 
(A, B) <r^ (C, D) ((£, E) «-> (G, //)) in each plaquette. The 
SU(2)-rotational and the exchange symmetries are emergent 
for the groundstate. For the SD phase in the spin-1/2 plaque¬ 
tte chain, as a result, the single groundstate has the emergent 
SU(2)-rotational and exchange symmetries. 


C. Modulated anti-ferromagnetic plaquette phase 

In contrast to the FP and SD phases, for the MAFP phase, 
a two-fold degenerate groundstate is detected. To discuss lo¬ 
cal magnetic properties of the states, let us first denote the 
two groundstate wavefunctions as with n e (1,2). For 
the system parameters 7' = 27 and h - 1.57, from the two 
degenerate groundstates, the local magnetizations 
are plotted at lattice sites in Fig. ©a). For the lo¬ 
cal magnetizations have a two-plaquette (eight-site) periodic 
structure where the first two sites, i.e., A and B, have zero- 
magnetizations = 0, the second two sites, i.e., C and D, 
have an opposite value of magnetizations depending on 7'/7, 
i.e., = -{Sf} > 0, the third two sites, i.e., E and F, have 

an maximum magnetizations {S^ - 1/2, and the fourth two 
sites, i.e., G and H, have an opposite value of magnetizations 
depending on 7'/7, i.e., {S^) - -{S^} < 0. For other values 
of system parameters, the characteristic behaviors of the local 
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FIG. 6: (color online) (a) Local magnetization (5^) at the lattice site 
i for y' = 2J and h = 1.57 in the MAFP phase, (b) Pictorial rep¬ 
resentation of the ground states with spin configuration. Red arrows 
indicate fully polarized spins. Two spins connected by a thick green 
line represents a singlet state, which means that the local magneti¬ 
zations at the spin sites are zero. Two spins connected by a thinner 
brown line represents an entangled pair states, where blue arrows de¬ 
note magnetizations at the spin sites. Note that each plaquette has an 
anti-ferromagnetic configuration of local magnetizations. 


magnetizations are not changed and only the values of the lo¬ 
cal magnetizations at the sites C, D, G, and H are determined 
by jyj. In this MAFP phase, also, the local magnetizations 
has a very characteristic property, i.e, an anti-ferromagnetic 
configuration in each plaquette. Note that for one-plaquette 
(four-site) shift, the local magentizations from |i/r 2 ) are equal 
to ones from \(j/\). This implies that one groundstate be¬ 
comes the other groundstate under one-plaquette (four-site) 
shift transformation. 

To obtain an explicit form of groundstates, we have dis¬ 
cussed the detailed properties of two-site spin correlations in 
the Appendix lAl From the discussions for it is found 
that for the first two sites i - A and B, the properties of 
the two-point spin correlations with the local zero magne¬ 
tizations satisfy the dimerized conditions in Eqs. dfiali-dfi^ 
and the conditions for a singlet state discussed in the SD 
phase. The spin state for sites A and B can then be written 
as = (ITa) Hb) - liA> ITb» / V2- For the third two sites 
i = E and F with the local maximum magnetizations, the 
two-point spin correlations between E/F and any site j in the 
system satisfy the conditions in Eqs. (l4al) and (l4bl) . which 
imply that each of the two sites is in a fully polarized state, 
i.e., IiAbf) = IT£)ITf)- For the second two sites (/ = C and 
D) and the fourth two sites (/ = G and H), the properties of 
the two-point spin correlations satisfy the dimerized condi¬ 
tions in Eqs. (I6all-(l6cli. Flowever, (‘^c/c'^c/c+i^ -1/4 even 

though (‘^c/c‘^c/G+ 1 ^ “ -1/4, which implies that the spin 
state for the dimerized two-sites is not in a singlet state. How¬ 
ever, other local correlations and magentizations discussed in 


the Appendix |B] allow us to write the spin state of the two 
sites as a linear combinations of two possible spin states, 
i.e., ^ with here i = C and 

G, where a, and b, are numerical coefficients depending on 
jyj and satisfying |a,p + = 1. Thus, the |i/'c/r) has 

a similar form with the \4/cd)- Comparing with the local 
magnetizations and the two-point spin correlations in the sec¬ 
ond two sites and the fourth two sites, the relations between 
the coefficients are given as - |hcl^ = “(*^0 “ l^cP) and 
flcl^cl = (^G\bG\ with |fl,p -H = 1, which leads to ac - Ihcl 

and ac - \bc\- Eor the second two sites (f = C and D) and the 
fourth two sites (/ = G and H), the spin states can be written as 
I^cd) - a |Tci£>)-|^l licTo) and |i/fc/r) = 1^1 ITci/rZ-a licT/r), 
where a and b are numerical coefficients depending on jyj 
with |ap + |hp = 1. Eor \ip 2 ), we have found similar properties 
of two-point spin correlations. Compared with the charac¬ 
teristic properties of local magnetizations and two-point spin 
correlations from \^\), it is found that the properties of them 
from |i/ri) are equal to ones from |i/r 2 ) for one-plaquette shift, 
which implies |i/' 2 ) are equal to under one-plaquette shift 
operation. Consequently, the groundstates li/'i) and |i/f 2 ) are 
given as 

l<Al) = ]~~[ [T 2 i) \(p2i,2i+\) |T 2 i+l,M) |T 2 i+l,(/) |^ 2 !+l, 2 i+ 2 ) , (8a) 

i 

I<A2) = ]~~[ |T2;,«) \'[2i,d)\v2i,2i+\) \X2i+l} |(^2i+l,2i+2) , (8b) 

i 

where = a|Ti>i;,/) - \b\\iiAjj)’ Ivij) = \b\\li,rij.i) - 

a |i,-,rT j,;), and \xi) = (|Ti,«) |Ti,rf>) / V2, Note that 

the two groundstates \ifri) and |i/' 2 ) are orthogonal to each 
other, i.e., {i/ri|i/f 2 ) = 0. 

The two groundstates in Eqs. (l8al) and dSbl) show the 
characteristic properties of their local symmetries as follows: 
(i) singlet states \x) are SU(2)-rotational invariant, (ii) fully- 
polarized states are polarized along the z-direction and then 
they are U(l)-rotational invariant on the x-y plane, and (iii) 
two-spin states \(p) and \(p) are U(l)-rotational invariant on the 
x-y plane because the x- and y-components of the local mag¬ 
netizations are zero, i.e., {S x) - 0 - {Sy}. Hence, for the 
MAEP phase, the two degenerate groundstates in Eqs. dS^ 
and (l8bl i are globally t/(l)-rotational invariant. Since the 
system Hamiltonian possesses the same t/(l)-rotational sym¬ 
metry with the two groundstates, for the MAEP phase, no 
global rotational symmetry is broken. Both of the two degen¬ 
erate groundstates are two-plaquette translational invariant, 
while our system Hamiltonian is one-plaquette translational 
invariant. Consequently, for the MAEP phase, the plaquette- 
translational symmetry breaking occurs and result in the de¬ 
generate groundstates. It is shown that the two degenerate 
groundstates can be understood within the Landau’s spon¬ 
taneous symmetry breaking picture. However, each of the 
two degenerate groundstates have the local SU(2) symmetry 
that cannot be explained within the Hamiltonian symmetry. 
Hence, the groundstates has such an emergent symmetry in 
the MAEP phase. 
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FIG. 7: (color online) (a) Local magnetization (5^) at the lattice site 
i for J[ = -1.417 and h = 0.3 in the MFP phase, (b) Pictorial repre¬ 
sentation of the groundstates with spin configuration. In contrast to 
the MAFP phase, note that each plaquette has a ferromagnetic con¬ 
figuration of local magnetizations in the MFP phase. 


D. Modulated ferromagnetic-plaquette phase 

The MFP phase is described by two-fold degenerate 
groundstates and |i/r 2 ). Let us discuss first about local 
magnetizations. We plot the local magnetizations from the 
two degenerate groundstates for 7' = -1.417 and h - Q3J 
in Fig. |7] The local magnetizations from |i/ri) have a two- 
plaquette (eight-site) periodic structure where the first two 
sites, i.e., A and B, have zero-magnetizations {S^} - 0, the 
second two sites, i.e., C and D, have (5^) = < 0, the 

third two sites, i.e., E and F, have an maximum magnetiza¬ 
tions (5^) = 1 /2, and the fourth two sites, i.e., G and H, have 
> 0. For other values of system parameters, 
the characteristic behaviors of the local magnetizations are not 
changed and only the values of the local magnetizations at the 
sites C, D, G, and H are determined by 7'/7. Compared with 
the local magnetizations from \ifri) for the MAFP phase in Fig. 
|6] one can notice that at the sites D and G, the local magneti¬ 
zations in the MFP phase have the opposite sign. This implies 
that the MFP phase has a distinguishable property of local 
magnetizations from the MAFP phase, i.e., a ferromagnetic 
configuration in each plaquette. For \il/ 2 }, also, this difference 
can be easily conformed from the local magnetizations from 
the MAFP phase in Fig. |6] and the MFP phase in Fig. [T] 
Furthermore, note that for one-plaquette (four-site) shift, the 
local magentizations from \ifr 2 } are equal to ones from 
Similar to the MAFP phase, in the MFP phase, one ground- 
state becomes the other groundstate under one-plaquette shift 
transformation. 

Actually, from the calculations of two-point spin correla¬ 
tions for |i/ri) in this MFP phase, we have found that for 
the first two sites i - A and B with the local zero magne¬ 
tizations and the third two sites i = E and F with the lo¬ 


cal maximum magnetizations, the characteristic properties of 
the two-point spin correlations are the same with those in the 
MAFP phase. For the sites, their spin states have the forms; 

= (ITa) liB) - liA> ITb)) / V2 and = Its) ITf)- For 
the second two sites (i = C and D) and the fourth two sites 
(i = G and H), the properties of the two-point spin corre¬ 
lations have been found to be the same with those of the 
fourth two sites (/ = G and H) and the second two sites 
{i = C and D) in the MAFP phase, respectively. This means 
that in the MFP phase, the spin states for the second and the 
fourth two sites are given as \i//cD) = 1^1 ITcin) - « licTo) and 
\^ch) - a ITci//) - \b\ IIgTh), where a and b are numerical 
coefficients depending on 7'/7 with |ap -i- \b'^ - 1. Simi¬ 
lar to the MAFP phase, in this MFP phase, we find that the 
physical properties from |^!r 2 ) are equal to ones from |i/ri) for 
one-plaquette shift. This implies \\p2) are equal to un¬ 
der one-plaquette shift operation. Consequently, for the MFP 
phase, we obtain the groundstates \tl/\) and \^2) as 

= \X2i) |'^2i,2i+l) |T 2 (+ 1 ,«) |T 2 i+l,rf) |(C’ 2 ,+l, 2 i+ 2 ) , (9a) 

i 

I<A2) = niT 2i,«> \hi.d} |(f> 2i,2i+ 1 ) [T2i+1) 1 4’2i+ 1,2i+2 ) • (9b) 

i 

These two groundstates \\}/\) and |^!r 2 ) are orthogonal to each 
other, i.e., (i/^i 11 /^ 2 ) = 0. One can also find that one-plaquette 
translational operation on |(!'i/ 2 ) leads to |i/' 2 /i)- Basically, the 
two degenerate groundstates in Eqs. ( l9al i and ( l9b] i have the 
same characteristic local symmetries with those in Eqs. ( l8al ) 
and dSbl) in the MAEP phase. Also, both of the two de¬ 
generate groundstates are two-plaquette translational invari¬ 
ant, while our system Hamiltonian is one-plaquette transla¬ 
tional invariant. Thus, for the MEP phase, the two degen¬ 
erate groundstates in Eqs. (l9al l and (l9b] l can be understood 
within the Landau’s spontaneous symmetry breaking picture, 
i.e., one-plaquette translational symmetry breaking results in 
the reduced symmetry of them comparing to the Hamiltonian 
symmetry. Similar to the MAEP phase, also, each of the two 
degenerate groundstates have the emergent local SU(2) sym¬ 
metry that cannot be explained within the Hamiltonian sym¬ 
metry in the MEP phase. 

E. Staggered bond phase 

In the SB phase, there are the two degenerate groundstates 
\il/\) and 11 /^ 2 ). Prom the two degenerate groundstates, for 
7' = 0.87 and h - 1.57, we plot the local magnetizations 
at lattice sites in Pig. [Sja) and the two-spin correlations as 
a function of lattice distance r between site i and / - 1 - r in 
Pig. |9] Por \^\) (I 1 A 2 )) in Pig.[8la), one can easily notice that 
the local magnetizations have one-plaquette periodic structure 
where the first (second) two sites have zero-magnetizations 

= 0 and the second (first) two sites have maximum mag¬ 
netizations {S- 112. We have found that the local mag¬ 
netizations are not changed for other values of system pa¬ 
rameters. Por the first two sites i - A and B, Pigs. |9|a) 
and|9jb) show that the properties of two-point spin correla¬ 
tions for are summarized as (i) = -1/4 and 
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FIG. 8: (color online) (a) Local magnetization (5^) at the lattice site 
i for Ji = 0.87 and h = 1.57 in the SB phase, (b) Pictorial represen¬ 
tation of the groundstates with spin configuration in the SB phase. 




(SPU = 0 for r > 1, and (ii) = 0. Also, 

) = 0 for a + a' (not displayed) have been ob¬ 
served. Since the local magnetizations are zero at the sites 
i - A and B, the properties of the two-point spin correla¬ 
tions imply that the two sites are dimerized and satisfy the 
conditions for a singlet state discussed in the SD phase, i.e., 
= (ITa) Hb) - Ha) ITb» / V2- Also, for the second two 
sites i = C and D, Figs. |9jc) and |9td) show that the two- 
point spin correlations for |(/ri) have the characteristic proper¬ 
ties, i.e., (i) = 1/4 = for j e {C,D} and 

= 0 for / 4 C and D, and (ii) = 0. We 

have noticed numerically that {S J ) - 0 for a + a' (not 
displayed). Since the local magnetizations have their maxi¬ 
mum value at the sites i = C and D, the properties of the two- 
point spin correlations imply that each of the two sites satisfies 
the conditions for a fully polarized state discussed in the FP 
phase, i.e., |i/'cd) = ITc)ITd). Similar discussions can be made 
for the other groundstate, As one can notice easily, the 
first (second) two sites in |t/ri) have the same physical proper¬ 
ties with the second (first) two sites in \i// 2 }. Consequently, for 
the SB phase, the two degenerate groundstates can be written 
as 

|(A2) = 

These two groundstates and 11 /^ 2 ) are orthogonal to each 
other, i.e., {iff\\i// 2 ) = 0. 

For the SB phase, the explicit forms of the two degener¬ 
ate groundstates in Eqs. (II Pal i and (llObl i can be represented 
pictorially in Fig. |8lb). One can see that the two degenerate 
groundstates are one-plaquette translational invariant and also 




FIG. 9: (color online) Two-point spin correlations with i = 

A,B,C,D and a = x, y, z as a function of lattice distance r between 
two sites i and i + r for 7( = 0.87 and h - 1.57 in the SB phase. 


have a global U(l)-rotational symmetry on the x-y plane. Our 
Hamiltonian has the one-plaquette translational and the U(l)- 
rotational symmetries of the two groundstates. Within the 
Landau’s spontaneous symmetry breaking picture, the occur¬ 
rence of the doubly degenerate groundstates for the SB phase 
could not be understood because no responsible symmetry for 
the degenerate groundstates exist in the system Hamiltonian. 
One can also notice that under the vertical-to-horizontal site 
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FIG. 10: (color online) (a) Local magnetization (5^) at the lattice site 
i for 7' = 1.43/ and h = Oin the SAFP phase, (b) Pictorial represen¬ 
tation of the groundstates with spin configuration in the SAFP phase. 
Note that each plaquette has an anti-ferromagnetic configuration of 
local magnetizations in the SAFP phase. 


exchange transformation, one groundstate becomes the other 
groundstate but the system Hamiltonian is not invariant. This 
fact implies that if the two degenerate groundstates are orig¬ 
inated from a spontaneous symmetry breaking, the vertical- 
to-horizontal site exchange symmetry might play a signihcant 
role to understand an emergent symmetry that is responsible 
for the two groundstates in the SB phase. Similar to the MAFP 
and MFP phases, also, each of the two degenerate ground- 
states have the emergent local SU(2) symmetry that cannot be 
explained within the Hamiltonian symmetry in the SB phase. 


F. Staggered anti-ferromagnetic plaquette phase 

Two degenerate groundstates are detected in the SAFP 
phase, as denoted as |t/ri) and \i// 2 }. From the two groundstates, 
we plot the local magnetizations at the lattice site 

i in Fig. [iota). In the SAFP phase. Fig. Q^a) shows that for 
the local magnetizations has a two plaquette (eight-site) 
periodic structure where (i) the first two sites, i.e., A and B, 
have a minimum magnetization (S^.) - -1/2, (ii) the sec¬ 
ond two sites, i.e., C and D, have (5^) = > 0^ (iii) 

the third two sites, i.e., E and F, have a maximum magneti¬ 
zation = 112, and (iv) the fourth two sites, i.e., G and 
H, {Sy) - —{Sy) < 0. For other values of system parame¬ 
ters, the characteristic behaviors of the local magnetizations 
are not changed and only the values of the local magnetiza¬ 
tions at the sites C, D, G, and H are determined by In 
this SAFP phase, also, the local magnetizations has an anti¬ 
ferromagnetic configuration in each plaquette. Note that for 
one-plaquette (four-site) shift, the local magentizations from 
|i/f 2 ) are equal to ones from \xjj\). This implies that one ground- 


state becomes the other groundstate under one-plaquette shift 
transformation. 

Actually, from the calculations of two-point spin correla¬ 
tions for li/ri) in this SAFP phase, we have found that for the 
first two sites i = A and B with the minimum magnetizations 
(the third two sites i - E and F with the maximum magneti¬ 
zations), the two-point spin correlations between AjB (EfF) 
and any site j in the system satisfy the conditions in Eqs. (l4al i 
and (l4bl i. which imply that each of the two sites is in a fully 
polarized state, i.e., \ijfAB) = IXa/IXb) and \i//EF) = ITe/ITf)- 
For the second two sites i - C and D (the fourth two sites 
i = G and H), the properties of the two-point spin correla¬ 
tions have been found to be the same with those of the second 
two sites i = C and D (the fourth two site i - G and H) in 
the MAFP phase, respectively. This means that in the SAFP 
phase, the spin state for the second (fourth) two sites has a 
same form of the spin state for the second (fourth) two sites 
in the MAFP phase, i.e., {(//cd) = alTcXo) - 1^1 IXcTd) and 
ken) = \b\ ITcXh) - a |XcT/r), where a and b are numerical 
coefficients depending on /'// with |ap + \b\^ - 1. Simi¬ 
lar to the MAFP phase, we hnd that in this SAFP phase, the 
physical properties from ^ 2 ) are equal to ones from ki) un¬ 
der one-plaquette shift transformation. Consequently, for the 
SAFP phase, we obtain the groundstates ki) and k 2 ) as 

kl) = ]””[ k2i,r;2i'+l,/) \^2i+l,ud} k2i+l,r;2i'+2,/), (1 la) 

i 

k2) = ]””[ |H'2/,Hrf) \4>2i,r,2i+l,l) |•^J•2/+1,H^^) |i,02/+l,r;2i'+2,/),(l lb) 
i 

where j|j. 2 /,wfif) = |X 2 /,uX 2 /,£/) and |tt 2 /,M£/) = |T 2 /,MT 2 /,fi?/- 

In the SAFP phase, both of the two degenerate ground states 
in Eqs. (II lal and (II Ibl l are global I/(l)-rotationaI invari¬ 
ant. Since the system Hamiltonian has the same global ro¬ 
tational symmetry with the two groundstates, the global ro¬ 
tational symmetry has nothing to do with the occurrence of 
the two groundstates. However, both of the two degenerate 
ground states in Eqs. (II 1 al l and (II Ibl l for the SAEP phase 
can be transformed from one to the other under one-plaquette 
translational or spin-flip transformations. Then, the plaque¬ 
tte translational symmetry and the spin-flip symmetry break¬ 
ings seem to result in the two degenerate groundstates. How¬ 
ever, note that for h - 0, the system Hamiltonian has both the 
one-plaquette translational symmetry and the spin-flip sym¬ 
metry but for h + 0, it has only the one-plaquette transla¬ 
tional symmetry. Hence, for h = 0 in the SAFP phase, the 
two-fold degenerate groundstates can be understood from the 
breaking of both the plaquette translational symmetry and the 
spin-flip symmetry. For h + 0, since the system Hamiltonian 
does not have the spin-flip symmetry, more symmetry than 
the symmetry that the Hamiltonian possesses should be in¬ 
volved to be broken for the two degenerate groundstates. This 
situation cannot be understood within the Landau’s sponta¬ 
neous symmetry breaking picture. Although, in order to un¬ 
derstand why the two groundstates exist in the SAFP phase 
for h 0, one may possibly consider the spin-flip symme¬ 
try as an emergent symmetry which can be broken together 
with a plaquette-translational symmetry, still a question why 
such an emergent symmetry should be broken together with 
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FIG. 11: (color online) (a) Local magnetization (5 j) at the lattice site 
i for J’. = — \AAJ and A = 0 in the SFP phase, (b) Pictorial represen¬ 
tation of the groundstates with spin configuration in the SFP phase. 
In contrast to the SAFP phase, note that each plaquette has a ferro¬ 
magnetic configuration of local magnetizations in the SFP phase. 


the plaquette-translational symmetry breaking must be left to 
be answered. 


G. Staggered ferromagnetic plaquette phase 

Similar to the SAFP phase, there are two degenerate 
groundstates detected in the SFP phase. The wavefunctions 
can be denoted as and \ip 2 )- From the two groundstates, 
we plot the local magnetizations at the lattice site 

i in Fig. fTTl a). In the SFP phase. Fig. fTTl al shows that for 
\xjj\), the local magnetizations has a two-plaquette (eight-site) 
periodic structure where (i) the first two sites, i.e., A and B, 
have a minimum magnetization {S - -1/2, (ii) the second 
two sites, i.e., C and D, have (5^) = < 0, (iii) the 

third two sites, i.e., E and F, have a maximum magnetiza¬ 
tion = 1/2, and (iv) the fourth two sites, i.e., G and H, 
(5p) = > 0- For other values of system parameters, 

the characteristic behaviors of the local magnetizations are not 
changed and only the values of the local magnetizations at the 
sites C, D, G, and H are determined by J[IJ. Comparing to 
the SAFP phase, this SFP phase shows a ferromagnetic con¬ 
figuration of the local magnetizations in each plaquette. Also, 
for one-plaquette shift, the local magentizations from \ip 2 ) are 
equal to ones from \il/\). Then, under one-plaquette transfor¬ 
mation, one groundstate becomes the other groundstate under 
one-plaquette shift transformation. 

Similar to the case of the SAFP phase, from the calculations 
of two-point spin correlations for \ip\) in this SFP phase, we 
have found that for the first two sites i - A and B with the 
minimum magnetizations (the third two sites i = E and F with 
the maximum magnetizations), the two-point spin correlations 
between A/B (EfF) and any site j in the system satisfy the 


conditions in Eqs. (l4al i and (l4ht . which imply that each of the 
two sites is in a fully polarized state, i.e., \i//AB) = IXa) IXb) and 
\^ef) - ITe/ITf)- However, for the second two sites i - C 
and D (the fourth two sites i - G and H), the properties of the 
two-point spin correlations have been found to be the same 
with those of the fourth two site i - G and H (the second 
two sites i = C and D) in the SAFP phase, respectively. This 
means that in this SFP phase, the spin state for the second 
(fourth) two sites has a same form of the spin state for the 
fourth (second) two sites in the SAFP phase, i.e., \iJ/cd) = 
1^1 ITcXD)-a IXcTd)j and Itf/cn) - a |TGX/r)-|^l IXcT/r) where a 
and b are numerical coefficients depending on 7'/7 with |ap -i- 
|hp = 1. Similar to the SAFP phase, we find that in this SFP 
phase, the physical properties from \^ 2 ) are equal to ones from 
under one-plaquette shift transformation. Consequently, 
for the SFP phase, we obtain the groundstates \^\) and \[f/ 2 ) as 

l<Al) = j~~[ l(£’ 2 i>; 2 i'+l,/) | 1 T 2 (+I,»rf) l 02 i+l,r; 2 i'+ 2 ,/),( 12 a) 

i 

I<A 2 ) = ]””[ |H' 2 /,H(/) l<^ 2 i,r; 2 )+l,;) \^2i+\,r,2i+2,l),i^'2^) 

i 

The groundstates in Eqs. ( I12al l and (I12bl ) in the SEP phase 
has a very similar form with the groundstates in Eqs. (II lal) 
and (II Ibl) in the SAEP phase. Then, the groundstates in the 
two phases have a very similar symmetry each other. Eor this 
SEP phase, one groundstate becomes the other groundstate 
for one-plaquette translation or spin-flip operation, which im¬ 
plies that the plaquette translational and the spin-flip symme¬ 
try breakings are responsible for the two degenerate ground- 
states. However, similar to the SAPP phase, the Hamilto¬ 
nian does not have the spin-flip symmetry h + 0 and then 
more symmetry than the symmetry that the Hamiltonian pos¬ 
sesses should be involved to be broken for the two degenerate 
groundstates. Thus, the occurrence of the two-fold degenerate 
groundstates in the SEP phase is not fully understood within 
the Landau’s spontaneous symmetry breaking picture. 

H. Anti-ferromagnetic plaquette phase 

The APP phase also have a two-fold degenerate ground- 
states. Prom the two groundstates, we plot the local magneti¬ 
zations {t^n\S at the lattice site i in Pig.[T2ja). The prop¬ 
erties of the local magnetizations in Pig. \TTl a) can be sum¬ 
marized as follows. Eor \^i) (li/fi)), the local magnetizations 
has a one-plaquette (four-site) periodic structure where (i) the 
first two sites, i.e., A and B, have a maximum (minimum) 
magnetization {S ^ = 1/2 (-1/2), (ii) the second two sites, 
i.e., C and D, have a minimum (maximum) magnetization 
{S^) - -112 (1/2). Eor other values of system parameters, 
the characteristic behaviors of the local magnetizations are not 
changed. Comparing to the local magnetizations in plaquettes, 
each plaquette has an anti-ferromagnetic configuration of the 
local magnetizations. Also, for one-plaquette shift, the local 
magentizations from |i/' 2 ) are equal to ones from \il/\). Then, 
under one-plaquette transformation, one groundstate becomes 
the other groundstate under one-plaquette shift transforma¬ 
tion. Prom the calculations of two-point spin correlations for 
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FIG. 12: (color online) (a) Local magnetization (5z) at the lattice site 
i for y' = -2|y| with / < 0 and = 0 in the AFP phase, (b) Pictorial 
representation of the groundstates with spin configuration in the AFP 
phase. Note that each plaquette has a ferromagnetic configuration of 
local magnetizations in the AFP phase. 


FIG. 13: (color online) (a) Local magnetization (S^) at the lattice 
site I for 7' = 2|y| with 7 < 0 and /; = 0 in the F phase, (b) Pictorial 
representation of the groundstates with spin configuration in the F 
phase. 


the groundstates, we have confirmed that for the first two sites 
/ = A and B with the maximum (minimum) magnetizations 
(the third two sites i = C and D with the minimum (maximum) 
magnetizations), the two-point spin correlations between A/Z? 
{CjD) and any site j in the system satisfy the conditions in 
Eqs. (l4al i and (l4bl) . which imply that each of the two sites is 
in a fully polarized state. Consequently, the two groundstate 
for the AFP phase can be respectively written as 

l<Ai) = (13a) 

i 

i<a2) = 

i 

Figure[T2jb) shows the two groundstates in the pictorial repre¬ 
sentation for the AFP phase. By comparing with the symme¬ 
try of the Hamiltonian, one can notice emergent symmetries 
for each groundstate. For instance, for a lattice rotation of 
each plaquette, each groundstate is invariant but the Hamilto¬ 
nian is not. Also, for a vertical-to-horizontal site exchange in 
each plaquette, e.g, (A, B) (C, D), the groundstate is invari¬ 
ant but the Hamiltonian is not. Hence, each groundstate has 
such emergent symmetries not belonging to the Hamiltonian 
symmetry. 

Furthermore, note that the two degenerate groundstate 
are transformed from one to the other under the spin-flip, 
the plaquette-rotational (one-site in each plaquette), or the 
vertical-to-horizontal site-exchange transformations. How¬ 
ever, for h = 0, the system Hamiltonian has only the spin-flip 
symmetry. Fven for h + Q, the system Hamiltonian does not 
have the spin-flip, the plaquette-rotational (one-site in each 
plaquette), and the vertical-to-horizontal site-exchange sym¬ 
metries. For = 0 in the AFP phase, hence, although the 
two-fold degenerate groundstates can be understood from the 


breaking of the spin-flip symmetry, more symmetries such as 
the plaquette-rotational (one-site in each plaquette) and the 
vertical-to-horizontal site-exchange symmetries should be in¬ 
volved to be broken for the two degenerate groundstates. For 
h + 0, since the system Hamiltonian does not have the spin- 
flip symmetry, no responsible symmetry for the two degen¬ 
erate groundstates exist in the system Hamiltonian. This is 
very similar case to the SB phase. Consequently, the two de¬ 
generate groundstates in the AFP phase cannot be understood 
within the Landau’s spontaneous symmetry breaking picture. 
Although, in order to understand why the two groundstates 
exist in the AFP phase for h ^ 0, one may possibly con¬ 
sider emergent symmetries such as the spin-flip, the plaquette- 
rotational (one-site in each plaquette), and the vertical-to- 
horizontal site-exchange symmetries that should be broken 
together, still the question why such emergent symmetries 
should be broken together must be an unanswerable. 


I. Ferromagnetic phase 

In the F phase, we detect two degenerate groundstates. 
From the two groundstates, we plot the local magnetizations 
{ifr„\S at the lattice site i in Fig. [T3l a). The local magneti¬ 
zations have {S^} = 1/2 from |i/fi) and (5j^) = -1/2 from 
Similar to the case of the FP phase, it has been observed that 
any site j in the system satisfies the conditions in Eqs. (l4al l 
and (l4bl i in the two-point spin correlations. Then, the state 
for the system is in a product state of the spin states of each 
site. According to the values of local magnetizations, the two 
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groundstates for the F phase can be respectively written as 

Wi) = (14a) 

i 

\<Pi) - j~| ■ (14b) 

i 

These two groundstates are presented pictorially for the F 
phase in Fig. fOl bl. 

Note that this F phase exist for h — 0. One can also easily 
notice that under a spin-flip transformation, one groundstate 
becomes the other groundstate and the system Hamiltonian 
is invariant. Hence, the two degenerate groundstates can be 
understood by the spin-flip symmetry breaking. However, by 
comparing with the symmetry of the Hamiltonian, one can no¬ 
tice emergent symmetries for each groundstate. Similar to the 
FP phase, each groundstate has more symmetries such as the 
lattice-rotation and the exchange symmetries. Consequently, 
each groundstate in the F phase has such emergent symmetries 
not belonging to the Hamiltonian symmetry. 

V. DISCUSSIONS ON EMERGENT SYMMETRY AND 
SPONTANEOUS SYMMETRY BREAKING 

In the previous section, for each phase, we have discussed 
how the explicit forms of groundstates can be extracted from 
the characteristic properties of the local magnetizations and 
the two-point spin correlations in our model. Emergent sym¬ 
metries of groundstates have been discussed. In this section, 
as a summary, we will discuss how degenerate groundstates 
can be related to spontaneous symmetry breaking in associa¬ 
tion with the Landau theory. What our results can suggest to 
understand beyond the Landau’s symmetry breaking mecha¬ 
nism will be also discussed. 

As we discussed in the introduction, if one assumes that 
two degenerate groundstates \^\) and |i/r 2 ), i.e., \^\) + |t/r 2 ), 
are obtained from any method, i.e., numerical or analytical 
calculations for a system Hamiltonian H, they satisfy H\{li„) - 
EgsVI^n) and can have a unitary transformation U connecting 
each other, i.e., - U |i/r 2 )- If the Hamiltonian is invariant 

under the unitary transformation, i.e., UHU^ - H and the uni¬ 
tary transformation U is related to an element of a subgroup 
g of the Hamiltonian symmetry group G, the two degenerate 
groundstates are originated from a breaking of a symmetry 
consisting of the subgroup g because of |^i) |i/^ 2 )- This 

case corresponds to the spontaneous symmetry breaking in 
the Landau theory. According to the spontaneous symmetry 
breaking mechanism in the Landau theory, then, the actual key 
symmetries associated with the two degenerate groundstates 
are presented in Lig. [T4]for the spin-1/2 plaquette chain, i.e., 
one-plaquette translational and spin-flip symmetries for h - 0 
and one-plaquette translational symmetry for h + 0. As we 
discussed in Sec. IV, the spontaneous breaking of spin-flip 
symmetry is shown to induce the two degenerate groundsates 
in the L phase and the spontaneous breaking of one-plaquette 
translational symmetry induces the two degenerate ground- 
states in the MALP and MLP phases. In the cases of SAPP 
(h - 0) and SPP ih - 0) phases, the spontaneous breaking 



FIG. 14: (color online) Relevant symmetries to degenerate ground- 
states in the Hamiltonian of Eq. l[T}. Gj, Gp, and G;, denote the 
symmetry groups of the dimer, plaquette, and field Hamiltonians, 
respectively. Then, the symmetry group of the Hamiltonian in Eq. 
O is G = Gjf]Gpf] Gh for h 0. For h = 0, it becomes 
G = GrfFlGp. V-to-H , T\-pjaquene^ ^2 Spin-flip , tiwd Rpjacjitene in¬ 
dicate the vertical-to-horizontal site-exchange symmetry, the one- 
plaquette translational symmetry, the spin-flip symmetry, and the 
plaquette-rotational symmetry, respectively. 


of both spin-flip and one-plaquette translational symmetries is 
applicable well to understand the two degenerate groundsates. 

However, if the Hamiltonian is not invariant under the uni¬ 
tary transformation U, i.e., UHU^ + H, the unitary transfor¬ 
mation U is not related to any element of any subgroup g of 
the Hamiltonian symmetry group G. Then a symmetry de¬ 
scribed by the unitary transformation does not belong to the 
Hamiltonian symmetry. In fact, such a situation has been ob¬ 
served in the SALP (h + 0), SLP (h + 0), AFP {h = 0), 
SB, and AFP (h + 0) phases in the spin-1/2 plaquette chain 
system. Based on the relations between the unitary transfor¬ 
mations and the Hamiltonian symmetry, to be more precise, 
we can categorize as for a given system parameter, (i) of all 
unitary transformations connecting degenerate groundstates, 
some are relevant to Hamiltonian symmetry and some are not 
relevant to Hamiltonian symmetry, and (ii) all unitary trans¬ 
formations connecting degenerate groundstates are not rele¬ 
vant to Hamiltonian symmetry. 

The first case (i) corresponds to the SALP {h + 0), SLP 
(h + 0), and AFP (h = 0) phases. In the SALP (h + 0) and 
SLP {h + 0) phases, that is, the two degenerate groundstates 
can transform from one groundstate to the other groundstate 
under the one-plaquette translational or the spin-flip trans¬ 
formations, but the Hamiltonian is invariant only under the 
one-plaquette translational transformation. Also, in the AFP 
{h - 0) phase, the two degenerate groundstates can transform 
from one groundstate to the other groundstate under the spin- 
flip, the plaquette-rotational, or the vertical-to-horizonal site- 
exchange transformations, but the Hamiltonian is invariant 
only under the spin-flip transformation. Then, these examples 
show that more symmetries than Hamiltonian symmetry can 
be broken for a spontaneous symmetry breaking. Such situa¬ 
tions cannot be explained fully by the spontaneous symmetry 
breaking mechanism in the Landau theory. For a complete ex¬ 
planation of occurring such degenerate groundstates in a view 
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of spontaneous symmetry breaking, if such more symmetries 
than Hamiltonian symmetry are considered as emergent sym¬ 
metries, they can be called an entailed-emergent symmetry 
that is broken together with a broken-symmetry connecting 
two degenerate groundstates within the Hamiltonian symme¬ 
try. 

The second case (ii) corresponds to the SB and AFP (h + 0) 
phases in the spin-1/2 plaquette chain system. In the AFP 
(h + 0) phase, that is, the two degenerate groundstates can 
transform from one groundstate to the other groundstate un¬ 
der the spin-flip, the plaquette-rotational, or the vertical-to- 
horizonal site-exchange transformations, but the Hamiltonian 
is not invariant. In contrast to the first case (i), moreover, 
no symmetry in the Hamiltonian symmetry is responsible for 
the two degenerate groundstates. In the SB phase, also, the 
vertical-to-horizonal site-exchange transformation connects 
one groundstate to the other groundstate, but the Hamiltonian 
is not invariant under the transformation and the Hamiltonian 
symmetry does not have any symmetry responsible for the two 
degenerate groundstates within the spontaneous symmetry 
breaking mechanism in the Landau theory. In order to explain 
such situations for a view of spontaneous symmetry breaking, 
if one can introduce an emergent symmetry that is respon¬ 
sible for degenerate groundstates, those cases can be called 
spontaneous emergent symmetry breaking. In this sense, such 
emergent symmetries can be called a to-be-broken emergent 
symmetry in order to distinguish from emergent symmetries 
that occur in degenerate groundstates, because for instance, 
each of the two degenerate groundstates in the F phase has 
emergent symmetries such as the vertical-to-horizontal site- 
exchange symmetry and the plaquette-rotational symmetry. 

However, although introducing such emergent symmetries 
associated with degenerate groundstates could explain occur¬ 
ring such degenerate groundstates in an extended view of the 
spontaneous symmetry breaking mechanism, the raised ques¬ 
tion, i.e., how such emergent symmetries associated with de¬ 
generate groundstates in a view of spontaneous symmetry 
breaking are related to Hamiltonian symmetry, is still left to 
be answered. To answer on the question, let us consider the 
Hamiltonian symmetry group in the spin-1/2 plaquette chain 
model in Eq. dB. In Fig. m we draw a schematic dia¬ 
gram to show the key symmetries responsible for degenerate 
groundstates in the symmetry groups of the dimer, the plaque¬ 
tte, and the held Hamiltonians. Gd, Gp, and Gh indicate the 
symmetry groups of the dimer, plaquette, and held Hamilto¬ 
nians, respectively. The symmetry group of the Hamiltonian 
in Eq. ([Til can be presented as a common subgroup of the 
symmetry groups Gd, Gp and Gh, i.e., G - GdClGpf^Gh for 
h + 0 and G = GdClGp for h - 0. Eor the spontaneous 
symmetry breakings in the E, MAER, MEP, SAEP (h - 0) and 
SEP (h = 0) phases, straightforwardly, the broken-symmetries 
are the one-plaquette translational or/and the spin-flip symme¬ 
tries belonging to the Hamiltonian symmetry group G. One 
may also easily notice that all so-called emergent symme¬ 
tries responsible for degenerate groundstates in the spin-1/2 
plaquette chain system belong to the common subgroups of 
pair of the symmetry groups Gd, Gp and Gh, i.e., G - G with 
G = (Grf n Gp) \JiGp n Gh) \J(Gh n Od). For instance, in the 


SB phase, the vertical-to-horizontal site-exchange symmetry 
belongs to the common group GdClGp - G. Furthermore, 
compared to the other phases, it should be noted that the two 
degenerate groundsates in the SB phase have a local SU(2) 
symmetry from the singlet states in Eqs. (llOab and (llObI) . Sim¬ 
ilarly, for in the MAFP (MFP) phase, a local SU(2) symmetry 
also appears in the two degenerate groundstates in Eqs. (l8al) 
and dSbl i (Eqs. ( l9al ) and (l9b]i). although the occurrence of the 
two degenerate groundstates can be explained by the spon¬ 
taneous breaking of the one-plaquette translational symmetry 
belonging to the Hamiltonian symmetry. Such an occurrence 
of the local SU(2) symmetry cannot be understood within the 
Hamiltonian symmetry group because only the dimer Hamil¬ 
tonian Hd has the local SU(2) symmetry. Hence, these facts 
imply that a symmetry belonging to the symmetry groups of 
the constituent Hamiltonians, i.e., Hd, Hp and Hh, plays a role 
for degenerate groundstates. In other words, for given system 
parameters, all symmetries belonging to the largest symmetry 
group Gd^jGpiJ Gh play a significant role for the system to 
reach its lowest energy states, i.e., degenerate groundstates. 
In a system, then, symmetry of degenerate groundstates might 
be determined within a largest common symmetry groups of 
constituent Hamiltonians, i.e., G, describing the system. Con¬ 
sequently, degenerate groundstates can be understood by an 
extended spontaneous symmetry breaking picture including 
spontaneous emergent symmetry breakings, and they do not 
have a broken-symmetry within a largest common symmetry 
of constituent Hamiltonians describing a given system but can 
have more symmetries than the largest common symmetry. 

VI. SPIN STRUCTURE FACTOR 

Experimentally, the quantum spin-1/2 plaquette chain 
can be realized in a number of magnetic compounds, es¬ 
pecially, some real insulating magnetic materials such as 
KCuL) 2 Dy][Mo(CN) 8 ] M, [Ee(H20)(L)][Nb(C^][Ee(L)] 
m and Dy(N 03 )(DMSO) 2 Cu(opba)(DMSO )2 lH. Those 
materials can be used to confirm an existence of emergent 
symmetries of groundstates. To do this, a way to detect the 
studies phases in our model system is to observe their char¬ 
acteristic magnetic properties, i.e., spin structure factors that 
can be observed by using neutron scattering experiments. Eor 
instance, for antiferromagnetic transverse-field Ising models 
in the pyrochlore lattice, a pinch point structure predicted the¬ 
oretically ll^[^ has been observed to correspond to a singu¬ 
larity in the spin structure factor of the spin-flip channel in the 
experiment of neutron scattering ll^ . In this section, thus, we 
will discuss spin structure factors for the eight phases in our 
Hamiltonian because the local magnetizations are zero in the 
SD phase. 

Eor one-dimensional lattice systems with N sites, a spin 
structure factor can be defined by the fourier transformation 
of spin correlation function. Eor z-direction, the spin structure 
factor can be defined as 

Z Z 

m=l j=l 
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where r = j - m is the lattice distance, q e [0,2n], and 
mji - I, - ■ ■ ,N. With the mapped one-dimensional chain 
structure in our system, we calculate spin structure factor Siq) 
to investigate the characterized behaviors of spin structures in 
the momentum space q. Actually, we have calculated S(q) 
numerically and analytically based on the two-point spin cor¬ 
relations studied in the previous sections for each phase. In 
Fig. [15] from numerical calculations, the spin structure fac¬ 
tor densities S{q) as a function of q with 0 < g' < 27r are 
plotted for four different phases with given parameters. It is 
shown that for each phase, the S{q) has a unique peak struc¬ 
ture for 0 < q < 2n\ (a) for the FP phase with the average 
magnetization = 1/2, there is a peak q — 0, (b) for the 
MAFP phase with — 1/8, seven peaks ?A q - knIA with 
k - 0,1,2,3,5,6 and 7, (c) for the SB phase with - 1/4, 
three peaks at ^ = knl2 with k -Q,l and 3, (d) for the SAFP 
phase with - 0, four peaks atq - knf A with k = 1,3,5 and 
7, (e) for the SFP phase with - 0, four peaks at q = knjA 
with k - 1,3,5 and 7, (f) for the MFP phase with \ /4, 
three peaks at q - knl2 with k - and 3, (g) for the AFP 
phase with - 0, two peaks at ^ = kn 12 with k - \ and 3, 
and (h) for the F phase with - +1/2, one peak at q - 0. 
Obviously, for the SD phase, Siq) - 0 because the local mag¬ 
netizations are zero for all lattice sites. Our analytic calcu¬ 
lation has given the same peak structure in the spin structure 
factor Siq) and is not presented. Such a characteristic peak lo¬ 
cation in Siq) then allows us to distinguish the different phases 
in our system. Consequently, the phases involving an emer¬ 
gent symmetry of groundstates can be observed by using a 
neutron scattering experiment. 


VII. CONCLUSIONS 

We have investigated a relation between degenerate ground- 
states and spontaneous symmetry breaking in the spin-1/2 
plaquette chain model. To do this, the iMPS representation 
with the iTEBD method are employed to calculate ground- 
states for a given parameter. The quantum fidelity has been 
shown to enable detecting degenerate groundstates for the sys¬ 
tem paramter space. By using the local magnetizations and the 
two-point spin correlations, the explicit forms of groundstates 
for each phase have been obtained. We have found that for 
whole parameter range, there are nine phases in which two 
phases (i.e., the FP and SD phases) have a single groundstate 
and seven phases have a two-fold degenerate groundstates. 

From unitary transformations connecting two degenerate 
groundstates each other, generally there are three types of de¬ 
generate groundstates: for a given system parameter, (i) all 
unitary transformations connecting degenerate groundstates 
are relevant to Hamiltonian symmetry, (ii) of all unitary trans¬ 
formations connecting degenerate groundstates, some are rel¬ 
evant to Hamiltonian symmetry and some are not relevant 
to Hamiltonian symmetry, and (iii) all unitary transforma¬ 
tions connecting degenerate groundstates are not relevant to 
Hamiltonian symmetry. The first case (i) corresponds to the 
F, MAFP, MFP, SAFP ih = 0), and SFP (/r = 0) phases, 
which is understood by the Landau’s symmetry breaking pic- 



q/ir 


qllT 


(C)' 


SB 

(d) 0-2 


SAFP 

1 1 






0.0 





0.0 0.5 1.0 1.5 2.0 

q/ir 


0.0 0.5 1.0 1.5 2.0 


(e) 


0.12 



0.0 0.5 1.0 1.5 2.0 

qllT 


0.0 0.5 1.0 1.5 2.0 

<ih 


(g) o-i 


0.0 0.5 1.0 1.5 2.0 

q/ir 



FIG. 15: (color online) Spin structure factors for (a) 7' = 0.27 with 
h = 1.57 in the FP phase, (b) 7/ = 27 with h = 1.57 in the MAFP 
phase, (c) J[ = 0.87 with h = 1.57 in the SB phase, (d) J[ - 1.437 
with h = Q m the SAFP phase, (e) 7/ = -1.447 with /r = 0 in the 
SFP phase, (f) 7,' = -1.417 with h = 0.37 in the MFP phase, (g) 
7/ = -2|7| with h = -|7| (7 < 0) in the AFP phase, and (h) 7/ = 2|7| 
with h - 0 iJ <0)m the F phase. 


ture. The SAFP ih + 0), SFP ih + 0), and AFP ih = 0) 
phases belong to the second case (ii), which cannot be under¬ 
stood fully by the Landau’s symmetry breaking picture. The 
third case (iii) are the SB and ALP ih + 0) phases, which 
is beyond the Landau’s symmetry breaking picture. Lurther- 
more, the groundstates can have more symmetries than the 
Hamiltonian symmetry. Lor instance, the appearances of the 
local SU(2) symmetry for the MALP, MLP, SD, and SB phases 
and of the emergent vertical-to-horizontal site-exchange and 
plaquette-rotational symmetries for the L phase are not ex¬ 
plained by the Landau’s symmetry breaking picture. Such rich 
emergent phenomena occurring beyond the Landau theory in 
our model suggest an equal footing theory as an extension of 
the Landau’s spontaneous symmetry breaking, i.e., degener¬ 
ate groundstates are induced by a spontaneous breaking of 
symmetries belonging to a largest common symmetry of con¬ 
tinent Hamiltonians describing a given system but can have 
more symmetries than the largest common symmetry. 

Linally, the characteristic properties of the spin structure 
factors in the different phases have been discussed. It is shown 
that the spin structures have the unique peak structures that 

























































16 


can be distinguished from one another. Such distinguishable 
peak structures can be observed by using a neutron scatter¬ 
ing experiment, which can be an experimental evidence for 
the extension of the Landau’s spontaneous symmetry break¬ 
ing theory. 
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Appendix A: Spin states for A,B,E and F sites in the MAFP 
phase 

In the Appendix, we will discuss about the two-point spin 
correlations in order to extract an explicit form of degen¬ 
erate groundstates in the MAFP phase. From the first two 
sites i - A and B in Figs. fTbl al and fTbl bl. the two-point 
spin correlations are shown that (i) = -1/4 and 

{S^SU = 0 for r > 1, (ii) - 0. Also, we 

have observed that = 0 for a + a' (not dis- 

played). For the third two sites i - E and F, in Figs. [T6lc) 
and fTbl' d). the two-point spin correlations are shown that (i) 
{S%^pS]) - {syxs)} and (ii) = 0. Also, we 

have observed that {S‘^ipS" ) = 0 for a + a' (not displayed). 
As we discussed in the main text, such properties of the two- 
point spin correlations lead to the explicit form of the spin 
states as \^ab) = (ITa) lia) - lU) ITb)) / V2 for the first two 
sites A and B, and liffEp) - ITfi) ITf) for the third two sites E 
and F. 


Appendix B: Spin states for C, D, G and H sites in the MAFP 
phase 

For the second two sites i - C and D in Fig. [17] the proper¬ 
ties of the two-point spin correlations for 7' = 2J and h - 1 .57 
are summarized as follows: (i) - -0.1118 and 

{S = 0 for r > 1. (ii) {S ^5 has an aperiodic struc¬ 
ture, i.e., except for the first period, {S ^5 has an eight-site 

periodicity. The only difference between the properties of the 
first period and the other periods is “ -1/4 and 

= -0-2 = <5^c><^^c.8.ri>withm = 1,2,.... 

(hi) = 0 for cr 4 a', (iv) = 0. (v) 

has an eight-site periodic structure. Also, we have 
observed that ^ 

= 0. From the summary, it should be noted 
that spin correlation, = (5^(5^^^) for r > 1 

and for r > 1. These proper¬ 

ties of the two-point spin correlations at C and D satisfy the 
dimerzation conditions in Eqs. dbal l. (l6bl i. and (jbc]). How¬ 
ever, + -1/4 even though = -1/4, 

which implies that the spin state for the two sites is not in 




r 




FIG. 16: (color online) Two-point spin correlations (5^5^^,.) as a 
function of lattice distance between two sites /? and /? -l- r for 7/ = 27 
and h= 1.57 in the MAFP phase. Here,/! 6 {A,B,E, F). 


a singlet state. Furthermore, for different system param¬ 
eters in the MAFP phase, the values of (5^), 

and change even though = -1/4 does not 

change. Also, from the numerical calculation, the values of 
(S^), and are observed to depend on only 

J'/J. Actually, for any parameter in the MAFP phase, the 
= -1/4 implies that the state for the two sites is a 
linear combination of two possible spin states, i.e., ITci/t) 
licTo). and then it can be written as lifrco) - a{jyj) |TciD) + 
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FIG. 17: (color online) Two-point spin correlations (5^/o‘^c/o+r) 
a function of lattice distance between two sites C/D and C/D + r for 
J[ = 2J and h = 1.57 in the MAFP phase. 


b{J'JJ)\ic'\D) with \a{J'JJ)\^ + \b{J'JJ)\^ = 1. In terms of 
the numerical coefficients a and b, the local magnetizations 
and the two-point spin correlation are expressed as (5^) = 

= (|fl(7'/7)P - \b{rjJt)l2 and = 

|a(7'/7)||h(7'/7)]|COS0/2 with a relative phase 6 between a 
and b, respectively, from the expression of \^cd)- Numeri¬ 
cally, the relations of (5^) = -{S\,) and (5^5^) = (5^5},) 
are manifested from Figs. |6]and (TT] Also, for all the system 
parameters in the MAFP phase, we have numerically found 
that {SISD = {S^^SD = -(1 + liSDY'W + 2<5^»'/2/4, 
which implies that the relative phase 9 does not depend on the 
system parameters and 9 - n. Then, one has a freedom to set 
that fl is a positive real number and b can be written as -\b\. 
Consequently, a best expression for the spin state of the two 
sites C andDcanbe |i/ccd) = a{J'JJ)\'\ciD)-\b(,J'JJ)\ licTz)>- 

For the fourth two sites i - G and H in Fig. [18] the proper¬ 
ties of the two-point spin correlations for 7' = 27 and /i = 1.57 
are summarized as follows: (i) - -0.1118 and 

{S = 0 for r > 1. (ii) {S ^5 has an aperiodic struc¬ 
ture, i.e., except for the hrst period, {S ^5 has an eight-site 

periodicity. The only difference between the properties of the 
first period and the other periods is = -1/4 and 

- -0-2 - With m = 1,2,.... 

(hi) = 0 for a ^ a', (iv) {Sfsfj = 0. (v) 

has an eight-site periodic structure. Also, we have 
observed that (vi) For a + a', 

(5 ^5 = 0. From the summary, it should be noted that 

spin correlation, = (‘^G^(‘^G+r) for r > 1 and 




FIG. 18: (color online) Two-point spin correlations as 

a function of lattice distance between two sites G/D and GjH -l- r for 
7/ = 27 and h = 1.57 in the MAFP phase. 


= (5 ^)(5 for r > 1. These properties of 
the two-point spin correlations at G and H satisfy the dimer¬ 
ization conditions in Eqs. (ihal i. ( l6b] i. and (16^ . However, 
+ -1/4 even though = -1/4, which 

implies that the spin state for the two sites is not in a sin¬ 
glet state. Furthermore, for different system parameters in 
the MAFP phase, the values of {S^}, and (5p^j) 

change even though - -1/4 does not change. 

Similar to the pair state of the two sites C and D, the in¬ 
variant = -1/4 of the two sites G and H for 

any parameter in the MAFP phase allows the state of the 
MAFP phase as a linear combination of two possible spin 
states, i.e., \i//ch} = a'(7'/7) ITgXh) + ^'(7'/7) |icT/r> with 
|a'(7'/7)p -H |h'(7'/7)p = 1, where the numerical coefficients 
have been conhrmed to be independent on the external mag¬ 
netic held. In terms of the numerical coefficients a' and b', 
the local magnetizations and the two-point spin correlation are 
expressed as = -<5^> - (|fl'(7'/7)|2 - |h'(7'/7)|2)/2 
and <5^5^) = (ScS^) = |fl'(7'/7)||h'(7z/-/)]|cos6»72 with 
a relative phase 9' between a' and b', respectively, from the 
expression of |i/'gh)- Numerically, the relations of (S^) = 
and (S^S^) - are manifested from Figs. |6| 

and [18] Also, for all the system parameters in the MAFP 
phase, we have numerically found that {S = {S ^5 = 

-(1 -I- 2(5g))^^^(l -H 2(5^))^^^/4, which implies that the rel¬ 
ative phase 9' does not depend on the system parameters and 
then 9’ = n. Then, one has a freedom to set that a' is a posi¬ 
tive real number and b' can be written as -|h'|. As a result, a 
best expression for the spin state of the two sites G and H can 
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be \if/GH) = a'iJ'JJ) ITci/f) - \b'{J'JJ)\ licTn)- Interestingly, 
the |i/fG/r) has a similar form with the \ipcD)- Comparing with 
the local magnetizations and the two-point spin correlations in 
the sites in Figs.|6] [T7]and [18] it should be noted that {S- 
(Sjj) = and<5J5J) = <5^5 J), which give the 


relations between the coefficients, i.e., -\b'^ - -(a'^ - |h'p) 

and a\b\ = a'\b'\. With + \b\^ = 1 and a'^ + |hf = 1, one 
can obtained the relations a' = \b\ and \b'\ - a. In terms of the 
coefficients a and b, the spin state of the two sites G and H can 
be expressed as Ii/'gh) = \b{J'JJ)\ ITgXh) - a(J'JJ) |iGT/r>- 
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